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Abstract 

Aloe vera (Asphodeloideae) is a medicinal plant in which useful secondary metabolites are plentiful. Among the 
representative secondary metabolites of Aloe vera are the anthraquinones including aloe emodin and chrysophanol, which 
are tricyclic aromatic quinones synthesized via a plant-specific type III polyketide biosynthesis pathway. However, it is not 
yet clear which cellular responses can induce the pathway, leading to production of tricyclic aromatic quinones. In this 
study, we examined the effect of endogenous elicitors on the type III polyketide biosynthesis pathway and identified the 
metabolic changes induced in elicitor-treated Aloe vera adventitious roots. Salicylic acid, methyl jasmonate, and ethephon 
were used to treat Aloe vera adventitious roots cultured on MS liquid media with 0.3 mg/L IBA for 35 days. Aloe emodin and 
chrysophanol were remarkably increased by the SA treatment, more than 10-1 1 and 5-13 fold as compared with untreated 
control, respectively. Ultra-performance liquid chromatography-electrospray ionization mass spectrometry analysis 
identified a total of 37 SA-induced compounds, including aloe emodin and chrysophanol, and 3 of the compounds were 
tentatively identified as tricyclic aromatic quinones. Transcript accumulation analysis of polyketide synthase genes and gas 
chromatography mass spectrometry showed that these secondary metabolic changes resulted from increased expression of 
octaketide synthase genes and decreases in malonyl-CoA, which is the precursor for the tricyclic aromatic quinone 
biosynthesis pathway. In addition, anti-inflammatory activity was enhanced in extracts of SA-treated adventitious roots. Our 
results suggest that SA has an important role in activation of the plant specific-type III polyketide biosynthetic pathway, and 
therefore that the efficacy of Aloe vera as medicinal agent can be improved through SA treatment. 
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Introduction 

Aloe vera (Asphodeloideae) is a medicinal plant in which useful 
secondary metabolites are abundant [1,2]. Anthraquinones, which 
represent one class of Aloe vera secondary metabolites, are tricyclic 
aromatic quinones. Among the naturally occurring anthraquinone 
derivatives, aloe emodin and chrysophanol are the major 
compounds [3]. The tricyclic aromatic quinones of aloe have 
been proposed to be synthesized via the type III polyketide 
biosynthesis pathway. Recently, novel plant-specific type III 
polyketide synthases (PKSs), octaketide synthase (OKS), PKS4, 
and PKS5 were isolated from Aloe arborescens, and their functions 
were examined in E. coli. The heterologously expressed enzymes 
produced SEK and SEK4b, which have an octaketide structure, 



from eight malonyl-CoAs, but SEK and SEK4b were found to be 
shunt products of the type II polyketide biosynthesis pathway and 
have not been detected in plants [4,5] (Figure 1). This suggested 
that these novel plant enzymes might potentially be associated 
with biosynthesis of natural tricyclic aromatic quinones in aloe, but 
it remains unclear whether these enzymes produce end products 
such as aloe emodin and chrysophanol in vivo. 

Identification of secondary metabolites and their derivatives in 
plants is a first step to determine the associated biosynthetic 
pathways; however, many remain to be identified because plants 
are especially rich in secondary metabolites with complex and 
varied structures [6] . Metabolomics approaches have ushered in a 
new era for elucidating the complicated secondary metabolism of 
plants [7,8,9,10]. Among these techniques, mass spectrometry can 
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Figure 1. Tentative biosynthesis mechanism for tricyclic 
aromatic quinones. (A) Production of SEK4 and SEK4b compounds 
catalyzed by OKS, PKS4, and PKS5 in £ coli. (B) The chemical structures 
of tricyclic aromatic quinone derivatives. 
doi:1 0.1 371 /journal.pone.0082479.g001 

contribute to deducing the structures of unknown compounds 
[11], and comprehensive profiling of phenolic compounds 
including anthraquinones has been carried out in Cassia and 
Rheum species using high-performance liquid chromatography- 
electrospray ionization-mass spectrometry (HPLC-ESI-MS) 
[12,13,14,15]. In addition, anthraquinone derivatives in rhubarb 
extract that were biotransformed by rat liver and intestinal 
bacteria were identified by liquid chromatography-electrospray 
ionization- tandem mass spectrometry (LC-ESI-MS/MS) [16,17]. 
In aloe, metabolite profiling was recently carried out for Aloe vera 
leaves at different developmental stages using gas chromatogra- 
phy-ion trap-mass spectrometry (GC-IT-MS) and ultra-perfor- 
mance liquid chromatography-quadruple time of flight-mass 
spectrometry (UPLC-QTOF-MS) [18]. 

In this study, we used GC-MS and UPLC-ESI-MS analysis to 
investigate the changes in primary and secondary metabolites 
induced by elicitor treatments of Aloe vera adventitious roots. Plant 
cell culture systems using elicitors have been effective alternatives 
for production of secondary metabolites, and many abiotic and 
biotic elicitors have been used to induce or enhance the 
biosynthesis of secondary metabolites by stimulating plant cellular 
stress responses [19,20]. Plant-originated signaling molecules such 



as salicylic acid (SA), methyl jasmonate (MJ), and ethylene as well 
as microbe-derived molecules such as polysaccharides, glycopro- 
teins, and inactivated enzymes have been utilized for elicitation of 
secondary metabolites [21,22]. 

Here, we tested the abilities of plant-derived elicitors to activate 
the type III polyketide biosynthesis pathway in order to improve 
the production of tricyclic aromatic quinones. We analyzed the 
changes in metabolic profile and anti-inflammatory activity in the 
extracts of elicitor-treated adventitious roots. This work shows that 
SA activates the type III polyketide biosynthesis pathway, resulting 
in improved production of tricyclic aromatic quinones and 
increased anti-inflammatory activity. It suggests that the plant 
specific-type III polyketide biosynthetic pathway is regulated by 
endogenous SA signaling, and that the efficacy of Aloe vera as a 
medicinal agent can be improved through SA treatment. 

Methods 

Chemicals and reagents 

All plant growth media and growth hormones were obtained 
from Duchefa (Haarlem, The Netherlands). Aloe emodin, 
chrysophanol, aloin, rhein, and emodin were purchased from 
Santa Cruz Biotechnology (California, USA). Malonyl-CoA, 
succinyl-CoA, acetyl-CoA, other chemicals, and solvents were 
from Sigma-Adrich (St. Louis, USA). All reference standards were 
dissolved in 100% methanol, and MJ and SA were dissolved in 
99.9% (v/v) ethanol. 2-chloroethylphosphonic acid (ethephon) was 
freshly dissolved in distilled water. Dissolved MJ, SA, and 
ethephon were sterilized through 0.45 |0,m membrane filters 
(Whatman, Tokyo, Japan) 

Optimization of suspension culture conditions and 
elicitor treatments 

Young shoots of Aloe vera were provided by Kim Jeong Moon 
Aloe Co. Ltd Jeju, Korea). Explants of young leaves were 
inoculated onto MS medium [23] supplemented with 30 g/L 
sucrose, 0.5 mg/L 1-naphthaleneacetic acid (NAA), and 0.7 g/L 
plant agar to induce adventitious roots as described previously 
[24]. After 3 weeks of culture, the adventitious roots were 
transferred into 100 mL Erlenmeyer flasks containing 30 mL 
MS liquid containing 30 g/L sucrose and 0.3 mg/L indole-3- 
butyric acid (IBA). The pH of medium was adjusted to 5.8 and 
then autoclaved at 121 °C for 15 min. The induced adventitious 
roots were incubated on an orbital shaker (50 rpm) at 25°C under 
constant light conditions (light intensity: 7 |J.E/m 2 s). Various 
concentrations of MJ, SA, and ethephon were added to the 
medium to treat 35-d-old adventitious roots. 

Estimate of adventitious root growth 

The adventitious roots were harvested every 7 d up to 42 d 
during growth. Fresh weight was determined to 0.05 g accuracy, 
and then dry weight (DW) was measured after lyophilization. The 
growth ratio was calculated as follows: DW of harvested roots - 
DW of inoculated adventitious roots (initial DW) /initial DW as 
previously reported [25]. 

Quantification of aloe emodin and chrysophanol in 
intracellular and extracellular 

The extraction of aloe emodin and chrysophanol from the 
culture medium was performed as previously reported with some 
modifications [26]. Harvested growth medium (30 mL) was 
supplemented with XAD-4 (0.07 g) and continuously agitated at 
125 rpm for 5 days at 25°C. Then, the XAD-4 was collected by 
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vacuum filtration and resuspended in 5 mL 100% ethanol 
followed by re-incubation for 5 days at 25°C with 125 rpm 
agitation. The extract in ethanol was concentrated to 500 |iL. The 
HPLC analysis conditions for aloe emodin and chrysophanol 
within the adventitious roots and culture media were described 
previously [24]. All samples obtained from the same condition 
were analyzed with three independent replicates and all experi- 
ments were repeated twice. 

Sample preparation for UPLC-ESI-MS 

For analysis of secondary and primary metabolites in cultured 
adventitious roots, lyophilized Aloe vera adventitious roots that were 
treated with 0, 500, 1000, and 2000 uM SA for 24 h (180 mg, 
n = 5) were extracted with 70% ethanol under sonication for 
30 min at room temperature. The extracts were centrifuged at 
14000xg for 10 min at room temperature and then completely 
dried under nitrogen gas. The dried residue was resuspended in 
200 uL 70% methanol and filtered through a 0.2 |J.m polytetra- 
fluoroethylene membrane filter (Whatman, Tokyo, Japan). Each 
sample was collected in quintuplicate from independent culture 
flasks and experiments were performed twice. 

Sample preparation for GC-MS 

The 50 uL ethanolic extracts were completely dried under 
nitrogen gas, and agitated with N-methyl-N-(trimethylsilyl) tri- 
fluoroacetamide +1% trimethylchlorosilane and pyridine (1:2) at 
60°C for 15 min. The solutions were then transferred into 2 mL 
glass vials coupled with micro-inserts (Agilent, Santa Clara, CA) 
and capped immediately. 

Conditions for UPLC-ESI-MS 

The extracts were analyzed with a Waters Alliance LC coupled 
to Quattro microMS on AQUITY UPLC BEH C 18 columns 
(2.1x100 mm, 1.7 urn). The UPLC-ESI-MS was carried out using 
mixtures of 0.1% aqueous formic acid (solvent A) and 100% 
acetonitrile (solvent B) at 0.1 mL/min flow rate and maintained at 
40°C. The elution program was as follows: 30% B at 0 min, 30% 
B at 5 min, 35% B at 10 min, 70% B at 35 min, 70% B at 45 min, 
and 100% B at 50 min. The injection volume was 10 |J,L and 
capillary voltages were adjusted to +4.0 kV for positive mode and 
to —4.0 kV for negative mode. MS/MS analysis was performed 
using collision energy from 5 to 60 eV in positive and negative 
modes. 

Conditions for GC-MS 

The primary metabolites in aloe extracts were analyzed using a 
6890 gas chromatograph (Agilent Technologies, CA, USA) 
coupled with a JMS-GC mate (Jeol, Tokyo, Japan). A DB-5 
column (30 mx0.25 mm I.D., 0.25 |J,m film thickness, HP) was 
used with helium (99.9999% He) as a carrier gas at a constant flow 
of 1 mL/min. The oven temperature was held at 60°C for 5 min, 
ramped to 320°C at a rate of 10°C/min and held for 10 min. One 
microliter sample was injected in split mode (10:1). The ionization 
energy was 70 eV in electron impact mode. The transfer line and 
ion source temperatures were set at 300°C and 300°C, respec- 
tively. After a 300 sec solvent delay, mass spectra were obtained at 
20 scans per second with a mass range of 55-600 m/ z. 

Data processing and multivariate analysis 

All raw data obtained from GC-MS were converted to ASCII 
format. The raw data were reduced into 6 sec buckets and 
normalized by the total sum of intensities as previously reported 
[27]. Raw data files obtained from UPLC-ESI-MS were exported 



to MZ-mine software version 2.1 and filtered through the 
Savitzky-Golay filter method to remove the noise. The baseline 
was then corrected and peaks were detected. Drifted retention 
time (RT) between replicated samples was adjusted through peak 
alignment. Finally RT was normalized to reduce the deviation of 
RT between peak lists, and data alignment was performed using 
RANSAC aligner. 

The peak lists resulting from GC-MS and UPLC-ESI-MS were 
evaluated using multivariate analysis with SIMCA-P 12.0 (Ume- 
trics, Umea, Sweden). Unsupervised principal component analysis 
(PCA) was performed, and supervised partial least square 
discriminant analysis (PLS-DA) and orthogonal partial least 
squares-discriminant analysis (OPLS-DA) were processed to 
compare each elicitor-treated condition and obtain differential 
metabolites from the elicitor conditions. The major metabolites 
that were differential between the elicitor conditions were 
regarded as variable importance in the project (VIP) list and m/ z 
having a cutoff score above the 1 was selected in the PLS-DA 
model. 

Statistical analysis of raw files obtained from GC-MS, UPLC- 
ESI-MS, and results of other experiments were carried out using 
Statistica Version 10 (StatSoft Inc., OK, USA), and the Tukey 
Honestiy Significant Difference and Least Significant Difference 
test were performed at probability level of 0.05. 

Metabolite identification 

The compounds were identified based on mass spectra and RT 
of authentic compounds. Peaks obtained from GC-MS were 
identified by comparison with spectra of the National Institute of 
Standards and Technology library on the basis of MS spectra. 
Among the m/z from UPLC-ESI/MS, major compounds 
obtained from multivariate or statistical analysis were tentatively 
identified by MS/MS spectra. 

Isolation and RNA expression analysis of OKS genes in 
Aloe vera 

Total RNA was extracted by a modified lithium chloride 
method according to previous reports [28]. The cDNA synthesis 
was conducted using the Maxime™ RT PreMix (Oligo (dT) 15 Kit 
(iNtRON, Sungnam, Korea) according to the manufacturer's 
protocol. Synthesized cDNA was diluted 1 /5 and used as template 
to isolate OKS genes and for real-time PGR. 

Full length OKS genes were isolated using primers designed from 
AaOKS (Accession number: AY567707) and AaPKS4 (Accession 
number: FJ536 166.1). PCR was performed as follows: 94°C for 
5 min and 35 cycles of 95°C for 30 s, 45°C for 30 s, and 72°C for 
1 min 30 s. The obtained PCR products were cloned into T-blunt 
cloning vector (Solgent, Daejeon, Korea) and sequenced using an 
ABI 3730 XL DNA Analyzer (Applied Biosystems, CA, USA) with 
Ml 3 Forward and reverse primers. The inserts of positive clones 
were amplified to obtain OKS candidate genes. 

To carry out real time PCR, the Ubiquitin gene isolated from Aloe 
vera was used to normalize the Ct values of target genes (Accession 
No. EF539181). Gene-specific primers were designed using the 
Primer 3 program [29] and specificity was confirmed by 
sequencing with an ABI 3730 XL DNA Analyzer. PCR reactions 
were carried out using a Light cycler 480 (Roche, Mannheim, 
Germany). The thermal cycling conditions were as follows: 95°C 
for 5 min and 40 cycles of 95°C for 15 s, 58°C for 10 s, and 72°C 
for 10 s. The primer sequences used to amplify each gene are 
summarized in Table SI. 

Gene expression levels were analyzed by reverse transcription 
quantitative PCR (RT-qPCR) using gene specific primer sets. 
Quantities of total RNA were normalized by comparison of band 
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Figure 2. Kinetic analysis of Aloe vera suspension culture. (A) 

Biomass accumulation of Aloe vera adventitious roots cultivated in MS 
liquid medium supplemented with 0.3 mg/L IBA. (B) Accumulation 
patterns of aloe emodin and chrysophanol in Aloe vera adventitious 
roots cultured on MS medium including 0.3 mg/L IBA. Each value is the 
mean of replicates and error bars indicate standard deviation. 
doi:1 0.1 371 /journal.pone.0082479.g002 

intensities for Ubiquitin. Thermal cycling conditions for RT-qPCR 
were as follows: 94°C for 5 min and 28 cycles of 94°C for 30 s, 
58°C for 30 s, and 72°C for 5 min. The amplified PCR products 
were separated on 2% agarose gel. 

Luciferase activity driven by COX-2, AP-1 and NF-kB 

The JB6 P+ mouse epidermal cell line which was kindly 
provided by Dr. Zigang Dong (University of Minnesota, Austin 
MN) [30] was cultured in monolayers on minimal essential 
medium (MEM) supplemented with 5% (v/v) fetal bovine serum 
(FBS) and 0.1% penicillin/streptomycin at 37°C in a humidified 
atmosphere of 5% C0 2 . TheJB6 P+ cells were stably transfected 
with COX-2, AP-1, or NF-kB luciferase reporter plasmid 
containing the G418 resistance gene and maintained in 5% 
FBS-MEM supplemented with 200 mg/mL G418. The cells were 
seed to 96-well plate, and the plates were incubated in a 5% C0 2 
incubator at 37°G. When cultured cells reached approximately 
80% to 90%confluency, the cells were starved with 0.1% FBS- 
MEM for 24 h. After that, cells were treated with various 
concentrations (0, 10, 20, 40, and 100 ug/mL) of Aloe vera 
adventitious root extracts (non-treated, 500 uM SA-treated, 
1000 uM SA-treated, and 2000 uM SA-treated) for 1 h followed 
by exposure to UVB (0.05 J/cm 2 ) and incubation for 4 h. JB6 P+ 
cells treated with UVB were disrupted with 100 |J,L lysis buffer 
[0.1 M potassium phosphate buffer (pH 7.8), 1% Triton X-100, 
1 mil DTT, and 2 mM EDTA]. The luciferase activity was 
measured using a luminometer (Luminoskan Ascent; Thermo 
Electron, Helsinki, Finland). 
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Figure 3. Effect of elicitor treatments on aloe emodin and 
chrysophanol production in Aloe vera adventitious roots. Aloe 
emodin and chrysophanol production in Aloe vera adventitious roots 
treated with SA (A), MJ (B), and ethephon (C) for 24 h. Data are 
represented as means of replicate samples ± standard deviation. 
Statistical analysis was carried out using the Tukey test (* p<0.05, 
** p<0.01). Asterisks indicate significant differences compared to aloe 
emodin and chrysophanol contents obtained from non-treated 
adventitious roots. 

doi:1 0.1 371/joumal.pone.0082479.g003 

Results 

Optimization of suspension culture conditions and 
elicitor effects on accumulation of aloe emodin and 
chrysophanol in Aloe vera adventitious roots 

We began by optimizing the suspension culture conditions for 
Aloe vera adventitious roots prior to elicitation. Three-week-old 
adventitious roots cultured from solid media were transferred into 
MS liquid media including 30 g/L sucrose, together with 0.5 mg/L 
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Figure 4. Time-course analysis of aloe emodin and chrysophanol production following elicitation. Content of aloe emodin (A, C, and E) 
and chrysophanol (B, D, and F) in adventitious roots and culture medium following non and elicitation with 2000 uM SA (A and B), 500 |iM MJ (C and 
D), or 500 \M ethephon (E and F). Data are represented as means of replicate samples ± standard deviation. Statistical analysis was carried out using 
the Tukey test (* p<0.05, ** p<0.01). Asterisks indicate significant differences compared to aloe emodin and chrysophanol contents obtained from 
adventitious roots before elicitation. 
doi:1 0.1 371 /journal.pone.0082479.g004 



IBA, indole-3-acetic acid, or NAA. After 4 weeks of culture, the 
adventitious roots were grew normally in MS media supplemented 
with IAA and IBA, but not in that supplemented with NAA (data 
not shown). We then tested the influence of IAA or IBA at 
concentrations of 0, 0.1, 0.3, 0.5, and 1.0 mg/L on adventitious 
roots, followed by examination of MS, 1/2 MS, 2MS, B5 [31], 
and SH [32] media (Table S2). After 35 d of cultivation, the 
maximum biomass of adventitious roots and maximum production 
of aloe emodin and chrysophanol, as representative secondary 
metabolites in Abe vera, were found in MS media supplemented 
with 0.3 mg/L IBA (Figure 2 and Table S2). 

Next, elicitation conditions were investigated using Aloe vera 
adventitious roots cultured under optimized suspension culture 
conditions. MJ, SA, and ethaphon were added in various 
concentrations to 35-d-old adventitious roots, and contents of 
aloe emodin and chrysophanol in the adventitious roots and in the 
media were measured, respectively. Maximum production of aloe 



emodin and chrysophanol occurred in adventitious roots treated 
with 1000-2000 uM SA (Figure 3). Time course analysis revealed 
that when adventitious roots were treated with various concen- 
trations of SA, MJ, and ethephon, accumulation of aloe emodin 
and chrysophanol in adventitious roots was increased by more 
than 10-11 and 5—13 fold at 24 h, respectively, and that in the 
growth medium rose at 24-72 h of 2000 uM SA treatment 
(Figure 4A and 4B). Treatment with 500 uM MJ and 500 uM 
ethephon also led to increases in aloe emodin and chrysophanol at 
24 h (Figure 3B and 3G). With 500 uM MJ, endogenous levels of 
aloe emodin and chrysophanol steadily increased and spiked at 
24 h with over 4-7 fold and 3-5 fold increases, respectively. 
Exogenous levels of these metabolites increased at 48-72 h 
following MJ treatment (Figure 4C and 4D). Treatment with 
500 uM ethephon induced endogenous aloe emodin and chryso- 
phanol at 12 h, with 5- and 4-fold increases, whereas it did not 
promote secretion of the these metabolites into the growth culture 
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Figure 5. Alteration of primary metabolites in Aloe vera adventitious roots in response to SA. (A) OPLS-DA of primary metabolites 
obtained from adventitious roots treated with 0, 500, 1000, or 2000 \\M SA analyzed by GC-MS. (B) Alteration of primary metabolites associated with 
TCA and glycolysis in response to SA elicitation. The level of malonyl-CoA decreased in a SA dose-dependent manner. Levels of other metabolites did 
not change in response to SA elicitation. G6P: Glucose-6-phosphate, F6P: Fructose 6-phosphate, PEP: phosphoenolpyruvate. 
doi:1 0.1 371 /journal.pone.0082479.g005 



medium (Figure 4E and 4F). SA treatment showed the most 
marked effect on production of aloe emodin and chrysophanol, 
and thus we investigated the responses of the tricyclic aromatic 
quinone biosynthesis pathway to SA elicitation. 

Primary metabolite analysis using GC-MS and LC-MS 

We speculated that SA might affect primary metabolites that 
play roles as intermediates in the biosynthetic pathways for 
secondary metabolites as well as induce considerable changes in 
secondary metabolites including aloe emodin and chrysophanol. 
To monitor alterations in primary metabolites including malonyl- 
CoA, which is a precursor for the biosynthesis of tricyclic aromatic 
quinones, extracts from Aloe vera adventitious roots treated with 0, 
500, 1000, and 2000 uM SA for 24 h were analyzed by GC-MS 
and LC-MS. Primary metabolites profiling analyzed by GC-MS 
did not reveal any difference upon treatment with different 
concentration of treated SA. The score plot of OPLS-DA showed 
that the SA-treatment replicates did not form separate clusters 
based on peak intensity (Figure 5A). Similar patterns were 




Figure 6. Effect of SA elicitation on transcript accumulation of 

OKS and OKSL-1. (A) Transcript accumulation of OKS and OKSL-1 at 6 h 
of 0, 500, 1000, and 2000 SA treatment relative to that of Ubiquitin. 
(B) Time course analysis of gene expression of OKS and OKSL-1 in the 
presence of 1000 |iM SA relative to that of Ubiquitin. Each value is the 
mean of replicates and error bars mean standard deviation. Statistical 
analysis was carried out using the Tukey test (* p<0.05, ** p<0.01). 
Asterisks indicate significant differences compared to control groups. 
doi:1 0.1 371 /journal.pone.0082479.g006 
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observed in the PCA and PLS-DA score plots (data not shown). 
Peaks associated with glycolysis and the tricarboxylic acid cycle 
(TCA), which is the principal pathway to form malonyl-CoA, were 
not significantly changed (Figure 5B and Table S3) and were not 
segregated clearly in the same multivariate analyses (data not 
shown). Interestingly, LC-MS analysis revealed that malonyl-CoA 
was remarkably decreased in SA-treated adventitious roots in a 
concentration-dependent manner, suggesting that consumption of 
malonyl-CoA is affected by biosynthesis of tricyclic aromatic 
quinones including aloe emodin and chrysophanol (Figure 5B and 
Table S2). Succinyl-CoA and acetyl-CoA, which are involved in 
the TCA cycle, could not be detected in either non-elicited or 
elicited samples. 

Isolation and gene expression of Aloe vera OKS genes 

The decrease in malonyl-CoA might result from activation of 
the type III polyketide biosynthesis pathway. To investigate 
changes in expression levels of the OKSs following SA treatment, 
the fuU length cDNAs for Aloe vera OKS (AvOKS) and OKS like-1 
(AvOKSL-1) were isolated from adventitious roots based on the 
sequences of Aloe arborescense OKS {AaOKS, Accession number: 
AY567707) and PKS4 (AaPKS4, Accession number: FJ536 166.1). 
Each AvOKS candidate gene had a 1212-bp open reading frame. 
Compared with AaOKS, AaPKS4, and AaPKS5 in Abe arborescense, 
the deduced amino acid sequences of AvOKS and AvOKSL-1 were 
90-99% identical and included conserved active sites such as the 
chalcone synthase (CHS) active sites (Met 147, Gly 221, Gly 226, 
and Pro 388), the catalytic triad of CHS (Cys 174, His 316, and 
Asn 349), gatekeepers (Phe 225 and Phe 275), and Gly 207, Leu 
266, and Val 351 (Figure SI and Table S4) [5]. 

Real-time PCR and RT-qPCR analysis showed that the 
expression of AvOKS and AvOKSL-1 was increased in proportion 
to the concentration of treated SA (Figure 6 and Figure S2). The 
transcripts for AvOKS and AvOKSL-1 were up-regulated more than 
6-fold at 6 h of 1000 |J.M SA elicitation. This suggests that SA 
induced the expression of AvOKS and AvOKSL-1, and in turn, the 
elevated enzyme activity accelerated the condensation of malonyl- 
CoA, resulting in increased production of tricyclic aromatic 
quinone derivatives including aloe emodin and chrysophanol. 

Global metabolite analysis 

We carried out UPLC-ESI/MS analysis in adventitious roots 
treated with 0, 500, 1000, or 2000 uM SA for 24 h to identify 
which metabolites along with aloe emodin and chrysophanol were 
induced by SA. Based on RT and mass-to-charge ratio (m/z), 
1850 and 634 peaks were obtained in positive and negative mode, 
respectively (Figure S3). According to multivariate results of the 
UPLC-ESI-MS data sets, the PCA and PLS-DA score plots 
showed clear segregation between elicitor-treated groups and the 
untreated group in positive and negative modes (Figure S4 and 
Table S5). The OPLS-DA analysis also revealed obvious 
differences among different SA-treated groups in positive and 
negative modes, indicating that metabolite alterations were 
influenced by SA concentration (Figure 7 and Table S5). The 
compounds in the UPLC-ESI-MS data sets induced by SA 
treatment were classified by one way analysis of variance (p<0.05) 
in comparison with the control to sort significantly increasing 
variables. We eliminated variables below 200 m/z because the 
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Figure 7. OPLS-DA score plot. OPLS-DA score plots of control (black), 
adventitious roots analyzed by UPLC-ESI-MS in positive (A) and in negath 
doi:1 0.1 371 /journal.pone.0082479.g007 

basic structure of tricyclic aromatic quinones has a molecular 
weight above 200 m/z. Consequently, we obtained 370 and 130 
variables, which were assigned using RT and m/ z, in positive and 
negative modes, respectively (Table S6). 

Identification of SA-induced metabolites using UPLC-ESI- 
MS/MS 

Significandy increased variables in positive and negative modes 
were analyzed by UPLC-ESI/MS/MS to obtain structural 
information for unknown compounds induced by SA elicitation. 
Tentative structures of 37 variables including chrysophanol (SA 
18) and aloe emodin (SA 19) were identified in positive and 
negative modes based on the Chemical Analysis Working Group 
(Table 1) [33]. These compounds mainly displayed losses of H z O 
(18 Da), CO (28 Da), CH 3 (15 Da), CH 3 COOH (60 Da), and 
COOH (45 Da), and in rare cases, degradation of CH 3 CN 
(41 Da) was detected (Table 1) [16,34]. Among the identified 
compounds, variables numbered SA 9, 12, 14, 18, and 28 were 
designated as VIPs that were capable of discriminating between 
the samples in PLS-DA models. 

Chrysophanol (SA 18) and aloe emodin (SA 19) were observed 
in adventitious roots treated with 1000~2000 uM SA. Chryso- 
phanol (SA 18) was detected at 36.4 min, predominately due to 
the elimination of a CO residue to produce 225 m/z. A MS/MS 
spectra of deprotonated aloe emodin (SA 19) was monitored at 
19.5 min and produced one fragment at 240 m/z by loss of CHO. 
These MS/MS fragment patterns of aloe emodin and chrysopha- 
nol were the same as those of the authentic compounds. 

A total of 3 compounds (SA 11, 15, and SA 16) were highly 
likely to be tricyclic aromatic quinone derivatives. The MS/MS 
fragment patterns of SA 1 1 and SA 1 5 were similar to that of 
emodin, and SA 16 showed a similar fragment pattern to that of 
aloe emodin. SA 1 1 appeared at 21.2 min and gave a 
deprotonated molecule at 269 m/z. Those MS/MS spectra 
produced daughter ions of 241 m/z [M-H-CO]", 213 m/z [M- 
H-2CO]", and 195 m/z [M-H-2CO-H 2 0] ". This fragment 
pattern was similar to that of emodin [13], but the RT was 
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different from that of the pure standard compound, implying that 
the compound was not emodin per se but might have a similar 
structure. The SA 15 compound eluted at 8.6 min and produced a 
deprotonated daughter molecule at 313 m/z. The base peak was 
further cleaved into 269, 241, 225, and 197 m/z ions. The 
269 m/z resulted from elimination of COOH adduct. The 
241 m/z ion was produced by disassociation of CO, followed by 
degradation of a hydroxyl group to form 225 m/z and a CO 
adduct to give 197 m/z. We tentatively determined that the 
COOH substitution might exist in a structure resembling emodin 
[13]. In addition, the spectrum of SA 16 showed a similar pattern 
to aloe emodin. The initial fragment arose from disassociation of 
CH 3 COCH 3 to produce 269 m/ z and cleavage of CHO to form 
240 m/z as the base peak. The corresponding structure is 
considered to resemble aloe emodin with a CH 3 COCH 3 
substitution. 

The anti-inflammatory activity of SA-treated Aloe vera 
adventitious roots 

Metabolite profiling analysis revealed that a number of 
metabolites including aloe emodin and chrysophanol were 
induced by SA elicitation. Tricyclic aromatic quinones such as 
aloe emodin and chrysophanol have been reported to possess anti- 
inflammatory activity [35] . Based on these results, we investigated 
whether SA treatment of adventitious roots led to enhancement of 
anti-inflammatory activity in UVB-treated mouse skin cells, which 
are a well-established cell line for screening anti-inflammatory 
agents. 

All extracts obtained from non-treated and SA-treated adven- 
titious roots maintained mouse skin cell viability at 10-100 |ig/ mL 
concentrations (Figure S6). Extracts obtained from 1000 |J,M SA- 
treated and 2000 uM SA-treated adventitious roots at 25-100 ug/ 
mL suppressed the UVB-induced promoter activity of COX-2 
(Figure 8). UVB-induced transactivation of MF-kB and AP I was 
strongly repressed by extracts from 500 uM SA-treated adventi- 
tious roots at 100 u\g/mL and 1000 uM SA-treated adventitious 
roots at 25 (Xg/mL, respectively (Figure 8). Although NF-kB was 
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Figure 8. Effect of extracts from SA-treated adventitious roots 
on UVB-induced expression in mouse skin cells. UVB-exposed 
JB6 P+ cells that were stably transfected with plasmids containing the 
luciferase reporter gene fused to the COX-2 promoter (A), the AP-1 gene 
(B), or the NF- kB gene (C) were incubated with extract from 0, 500, 
1 000, or 2000 |rM SA-treated adventitious roots for 1 h and harvested 
after 4 h. Data are represented as means of replicate samples ± 
standard deviation. Statistical analysis was carried out with the Tukey 
test (* p<0.05, ** p<0.01). Asterisks indicate significant differences 
compared to groups treated with UVB alone. Control indicates the 
extract from SA-untreated adventitious roots. 
doi:1 0.1 371 /journal.pone.0082479.g008 

slightly inhibited by the extract from non-treated adventitious 
roots at 100 ug/mL, the extract could not suppress promoter 
activity of COX-2 and transactivation of AP I. These results 
indicated that anti-inflammatory activity was enhanced in the 
adventitious root extracts by SA elicitation, resulting from the 
induction of a number of metabolites with anti-inflammatory 
activity. In addition, anti-inflammatory activity was dependent on 
the concentration of SA used to treat adventitious roots. Although 
we could identify the metabolites induced by SA treatment as 
described above, further detailed study is needed to confirm which 
metabolites account for the anti-inflammatory activity. 



Discussion 

In this study, we investigated metabolic changes and enhance- 
ment of anti-inflammatory activity in SA-treated Aloe vera 
adventitious roots. Prior to this series of experiments, there was 
a need for development of in vitro cell culture systems for Aloe vera in 
order to retain cell lines in a controlled environment. Although 
previous studies attempted to optimize callus induction on solid 
media in Aloe species [36,37], suspension culture had been limited, 
because phenolic compounds released from the cultured cells 
eventually led to cell death [38]. In this work, we optimized 
suspension culture conditions for Aloe vera adventitious roots and 
overcame heavy browning (Figure 2 and Table S2). Using the 
adventitious roots, we examined the elicitation effects of MJ, SA, 
and ethephon on the production of aloe emodin and chrysopha- 
nol. 

MJ, SA, and ethephon are plant-derived elicitors that mediate 
the signal transduction involved in plant defense responses [39]. 
SA is associated with pathogen-related defense mechanisms and is 
required for establishment of plant systemic acquired resistance. 
Treatment with exogenous SA was previously reported to induce 
secondary metabolites, e.g. anthraquinones in Rubia cordifolia [39], 
soluble phenolic compounds in Matricaria chamomilla and Salvia 
miltiorrhiza [40,41], podophyllotoxin in Linum albums [42] and 
artemisinin in Artemisia annua L [43] . In our study, the production 
of aloe emodin and chrysophanol was increased by all of these 
elicitors, but SA remarkably elevated the level of aloe emodin and 
chrysophanol by more than 10-11 and 5-1 3 fold at 24 h, implying 
defense-related function of the quinone compounds (Figure 3 and 
Figure 4). 

We performed UPLC-ESI/MS/MS analysis to identify addi- 
tional SA-induced tricyclic aromatic quinones along with aloe 
emodin and chrysophanol. We identified 37 novel compounds 
induced by SA elicitation through MS/MS fragmentation, and 
three compounds were determined likely to have tricyclic aromatic 
quinone structures (Table 1). Our analysis of overall primary 
metabolite profiles and transcripts of OKS genes provided 
conclusive information about key components in the type III 
polyketide biosynthesis pathway that was previously uncertain in 
plants. In previous reports, it was predicted that tricyclic aromatic 
quinones such as aloe emodin and chrysophanol were synthesized 
from the precursor of malonyl-CoA in a process mediated by 
OKSs [4,5] . However the activity of the plant OKS enzymes was 
tested in E. coli, and it remained unclear whether the condensation 
process could indeed produce the end products in aloe, as opposed 
to the shunt products SEK and SEK4b (Figure 1). Our results 
showed that aloe emodin and chrysophanol biosynthesis correlated 
with OKS gene expression and malonyl-CoA content, supporting 
the hypothesis that the condensation of malonyl-CoA mediated by 
OKSs results in production of tricyclic aromatic quinones in Aloe 
vera (Figure 3,4,5,6, and Table S3). On the other hand, some 
tricyclic aromatic quinones were not induced following SA 
elicitation. Rhein and emodin were decreased or not changed 
following SA elicitation, respectively (data not shown), and 
indicating that SA elicitation might also affect the downstream 
pathway of the condensation process. 

Our metabolite analysis suggests that there is cross-talk in SA 
and MJ signaling in Aloe vera. UPLC-ESI-MS based-analysis 
showed clear differences between 500 u,M SA-treated and 500 uM 
MJ-treated groups in positive mode (Figure S5 and Table S5). This 
reveals that SA and MJ can lead to different responses in overall 
metabolite accumulation in Aloe vera adventitious roots. In many 
cases, SA and MJ have been found to act antagonistically [44], but 
some cooperative activity was observed in naphthodianthrone and 
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phloroglucinol production in Hypericum species and ginsenoside 
accumulation in Panax ginseng [45,46]. In Aloe vera, the biosynthesis 
of many secondary metabolites was SA- or MJ-dependendy 
regulated, but some metabolites such as tricyclic aromatic 
quinones seemed to be cooperatively regulated. 

Finally, we examined whether SA elicitation enhanced the anti- 
inflammatory activity of adventitious root extracts. We verified 
anti-inflammatory effects using JB6 P+ cells and the UVB model, 
which is suitable for examining chemopreventive effects of 
phytochemicals [47,48,49]. Our results showed that anti-inflam- 
matory activity was affected differently by extracts from adventi- 
tious roots treated with different concentrations of SA (Figure 8). 
This indicates that SA-induced metabolites including aloe emodin 
and chrysophanol might regulate anti-inflammatory activity in 
UVB-exposed mice skin cells. Our metabolite profiling data 
revealed that the identified compounds mainly possessed hydroxy, 
hydroxymethyl, and carboxyl groups. There is accumulating 
evidence that the presence of these residues might have an effect 
on inflammation or angiogenesis associated with inflammation 
[50,51,52]. For example, the presence of 2 hydroxy groups located 
in aloe emodin might play an important role in its anti- 
inflammatory activity [51]. Emodin with hydroxyl group at C-6 
position, rhein with carboxyl groups at C-6 position, and aloe 
emodin with hydroxymethyl group at C-3 position might 
contribute those anti-angiogenetic properties [50,51]. These types 
of evidence might partially account for the enhancement of anti- 
inflammatory activity in SA-treated adventitious roots. 

Abe vera is a well-known commercial crop; however, only a few 
metabolomic approaches have been applied to it. In this study, 
based on an optimized cell culture system for Aloe vera, we showed 
that SA elicitation led to activation of the tricyclic aromatic 
quinone biosynthesis pathway and the accumulation of secondary 
metabolites. Moreover, the extracts from SA-treated Aloe vera had 
enhanced anti-inflammatory activity in UVB-treated mouse skin 
cells. Taken together, these results provide a possible biological 
function of quinones related to SA-dependent defense responses 
and reveal the potential of SA-induced metabolites as chemopre- 
ventive agents. 

Supporting Information 

Figure SI Alignment of amino acid sequences of Aloe 

vera OKS and OKSL-1. OKS and OKSL-1 sequences from 
Aloe vera were compared with OKS from Aloe arborescens (Accession: 
AY567707.1), PKS4 from Aloe arborescens (Accession: FJ536166.1), 
and PKS5 from Aloe arborescens (Accession: FJ536 167.1). 
(TIF) 

Figure S2 Expression levels of OKS and OKSL-1 in 
response to SA treatment. RT-qPCR expression profile of 
OKS and OKSL-1 at 6 h of 0, 500, 1000, and 2000 uM SA 
treatment (Lane: Cont, SA500, SA1000, and SA2000) and time 
course analysis of gene expression of OKS and OKSL-1 in the 
presence of 1000 uM SA (Lane: 0, 3hr, 6hr, 12hr, and 24hr). 
Quantities of total RNA were normalized by comparison with the 
band intensity for Ubiquitin, and the PGR products for Ubiquitin, 
OKS, and OKSL-1 were separated on 2% agarose gels. 
(TIF) 
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Figure S3 Chromatograms obtained from UPLC-ESI- 
MS in positive and negative modes. 

(TIF) 

Figure S4 PCA and PLS-DA score plots. The PCA (A and 

C) and PLS-DA (B and D) score plots of control (black), 500 (red), 
1000 (blue), and 2000 (green) |iM SA-treated adventitious roots 
analyzed by UPLC-ESI-MS in positive (A and B) and negative (C 
and D) mode. 
(TIF) 

Figure S5 PCA, PLS-DA, and OPLS-DA score plots. PCA 

(A and D), PLS-DA (B and E), and OPLS-DA (C and F) score 
plots of control (black), 500 uM SA (blue), and 500 uM MJ (red)- 
treated adventitious roots analyzed by UPLC-ESI-MS in positive 
(A, B, and C) and negative (D, E, and F) mode. 
(TIF) 

Figure S6 Effect of extracts obtained from elicitor- 
treated adventitious roots on JB6 P+ cell viability. 

Extracts were obtained from Aloe vera adventitious roots untreated 
(Control) or treated with 500 uM SA, 1000 uM SA, or 2000 \lM 
SA. JB6 P+ cells were treated with the indicated amounts of each 
extract for 4 h, and then 20 uL CellTiter 96 Aqueous One 
solution was added to the cells and they were incubated for an 
additional 4 h. Cell viability was subsequently measured at 492 
and 690 nm. Data are represented as means of replicate samples 
± standard deviation. Statistical analysis was carried out using the 
Tukey test (* p<0.05, ** p<0.01). Asterisks indicate significant 
differences compared to control groups. 
(TIF) 
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(DOCX) 

Table S2 Effect of plant hormones and media on growth 
of Aloe vera adventitious roots and accumulation of aloe 
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PKS4, and PKS5 from Aloe arborescens. 

(DOCX) 

Table S5 Statistical parameters of PCA, PLS-DA, and 
OPLS-DA in positive and negative modes. 
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